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We investigate the entanglement entropy in gravity duals of confining large gauge theories 
using the proposal of pQ. Dividing one of the directions of space into a line segment of length / 
and its complement, the entanglement entropy between the two subspaces is given by the classical 
action of the minimal bulk hypersurface which approaches the endpoints of the line segment at 
the boundary. We find that in confining backgrounds there are generally two such surfaces. One 
consists of two disconnected components localized at the endpoints of the line segment. The 
other contains a tube connecting the two components. The disconnected surface dominates the 
entropy for / above a certain critical value /crit while the connected one dominates below that 
value. The change of behavior at / = /crit is reminiscent of the finite temperature deconfinement 
transition: for / < /crit the entropy scales as N^, while for / > lent as N^. We argue that a similar 
transition should occur in any field theory with a Hagedorn spectrum of non-interacting bound 
states. The requirement that the entanglement entropy has a phase transition may be useful in 
constraining gravity duals of confining theories. 
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1 Introduction 

Consider a d + 1 dimensional quantum field theory (QFT) on IR'^'*'^ in its vacuum state |0). 
Divide the d dimensional space into two complementary regions, 

B = IR^-i X (IR - Ii) , (1) 

where /; is a line segment of length /. The entanglement entropy between the regions A and B 
is defined as the entropy seen by an observer in A who does not have access to the degrees of 
freedom in B, or vice versa (see e.g. [2] for a recent review and references to earlier work). It 
can be calculated by tracing the density matrix of the vacuum, po = |0)(0|, over the degrees of 
freedom in B and forming the reduced density matrix 

Pa = TtbPo ■ (2) 
The quantum entanglement entropy Sa is then given by 

Sa = -TiaPa'^'^Pa ■ (3) 

The above construction can be generalized in a number of ways. In particular, one can replace 
the vacuum state |0) by any other pure or mixed state, and choose the submanifold of IR'^, A, to 
be different than ([1]). In this paper we will restrict to the choices above, which are sufficient for 
our purposes. 

The entanglement entropy ([3]) is in general UV divergent. To leading order in the UV cutoff 
a it scales like [H H] 

where Vd-i is the volume of IR"^"^ in ([T]). Note that (jlj) is independent of /. This turns out to 
be a general feature - the entropy is defined up to an / independent (infinite) additive constant. 
In particular, OiSa and differences of entropies at different values of / approach a finite limit 
as a — >■ 0. In d + 1 dimensional CFT with d > 1, the finite /-dependent part of the entropy is 
negative and proportional to Vd-i/l'^~^, while for d = 1 it goes like In/. 
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If the QFT in question has a gravity dual [5], it is natural to ask whether the entanglement 
entropy can be calculated using the bulk description. This problem was addressed in [1]. For 
the case of d + 1 dimensional large Nc conformal field theories with AdSd+2 gravity duals, the 
authors of [1] proposed a simple geometric method for computing the entanglement entropy and 
subjected it to various tests. This method is to find the minimal area d-dimensional surface 7 
in AdSd+i such that the boundary of 7 coincides with the boundary of A, which in the case ([1]) 
consists of two copies of IR'^"^ a distance / apart. The quantum entanglement entropy between 
the regions A and B is proportional to the classical area of this surface, 



where G)^ is the d+2 dimensional Newton constant and G\^^j^ is the induced string frame metric 
on 7. Note that the surface 7 is defined at a fixed time and ([5]) gives the entanglement entropy 
at that time. For static states, such as the vacuum, the resulting entropy is time independent o 
Also, since 7 is extended in the transverse IR'^"^ ([1]), the entropy ([5]) is proportional to its volume 
Vd-i- Thus, in this case it is better to consider the entropy per unit transverse volume. 

In non-conformal theories, the volume of the 8 — d compact dimensions and the dilaton are 
in general not constant. A natural generalization of ([5]) to the corresponding ten dimensional 
geometries is [H [7] 



The entropy is obtained by minimizing the action ([6]) over all surfaces that approach the boundary 
of y4 ([1]) at the boundary of the bulk manifold, and are extended in the remaining spatial 
directions. Since = Svr^a'^^f^, this gives an answer of order A*"^ in the 't Hooft limit Nc ^ 00 
with QgNc held fixed. 

It was shown in [1] that for AdS^ the prescription ([5]) successfully reproduces the known form 
of the entanglement entropy in two-dimensional conformal field theory, and that it gives sensible 
results when apphed to some higher dimensional vacua, such as AdS^ x S^. Nevertheless, some 
aspects of the proposal are not well understood. In particular, it is not clear how to extend it 
beyond leading order in l/N^.. 

In this paper we apply the proposal of [H [7] to confining backgrounds, such as [10], [H] . One 
of the motivations for this investigation is to subject the proposal ([6]) to further tests. Another is 
to study the / dependence of the entanglement entropy, which is in general difficult to determine 
in strongly coupled field theories. 

Gravitational backgrounds dual to confining gauge theories typically have the following struc- 
ture. As one moves in the radial direction away from the boundary, an internal cycle smoothly 
contracts and approaches zero size at the infrared (IR) end of space. The radial direction together 
with the shrinking cycle make a type of cigar geometry, with the IR end of space corresponding 
to the tip of the cigar. 

We will see that in such geometries there are in general multiple local minima of the action 
([6]) for given /. One of those is a disconnected surface, which consists of two cigars extended 
in IR''"^ and separated in the remaining direction in IR'^ by the distance /. A second one is a 

^Generalizations of the proposal of [T] to time dependent states were discussed in [5]. 




(5) 




(6) 
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connected surface, in which the two cigars are connected by a tube whose width depends on /. 
Since the two geometries are related by a continuous deformation, there is a third extremum of 
the action between them, which is a saddle point of ([6]). A natural generalization of the proposal 
of [1] to this case is to identify the entanglement entropy with the absolute minimum of the 
action. We will see that this leads to a phase transition in the behavior of the entanglement 
entropy as a function of /. 

For the disconnected solution, Sa dS]) does not depend on /. As mentioned above, the actual 
value of Sa depends on the UV cutoff, but if we are only interested in differences of entropies, or 
the derivative of the entropy with respect to /, we can set it to zero. For the connected solution, 
Sa depends non-trivially on /. For small I, it is smaller than that of the disconnected one. Thus, 
it dominates the entropy IQ. For / > /crit the action of the connected solution becomes larger 
than that of the disconnected one, and it is the latter that governs the entropy. Thus, in going 
from / < /crit to / > /crit, diSA gocs from being of order A'^^ to being of order N^. One can think of 
this change of behavior as a phase transition which, as we show, is typical in large Nc confining 
theories. 

Similar transitions between connected and disconnected D-brane configurations play a role in 
other contexts. In [8] an analogous transition is responsible for screening of magnetic charges in 
confining gravitational backgrounds; in [9] it governs the pattern of metastable supersymmetry 
breaking vacua in a brane construction of supersymmetric QCD. An important difference is 
that in all these cases the transitions involve the rearrangement of real branes, whereas the 
hypersurface whose area is being minimized here does not seem to have such an interpretation. 

The plan of the rest of the paper is as follows. In section 2 we present a general analysis 
of a class of gravitational backgrounds that arises in the construction of holographic duals of 
confining gauge theories. We show that in this class there are multiple local minima of the 
action (jS]), as discussed above. With some mild assumptions, we also show that for small / the 
global minimum of the action corresponds to a connected solution, while for large / it corresponds 
to a disconnected one. We also show that the connected solution does not exist for sufficiently 
large /. 

In sections 3 - 5 we illustrate the discussion of section 2 with a few examples. Section 3 
contains an analysis of the geometry of Nc D4-branes wrapped around a circle with twisted 
boundary conditions for the fermions. For QgN^ <C 1 this system reduces at low energies to pure 
Yang-Mills (YM) theory, while for QgNc ^ 1 it can be analyzed using the near-horizon geometry 
of the D4- branes [10]. In section 4 we describe the analogous D3-brane system, which for 
QsNc <^ 1 gives rise to YM in 2 + 1 dimensions. Some of the results of this section were obtained 
already in [7]. Section 5 contains an analysis of the warped deformed conifold (KS) background 
[TT] . which corresponds to a cascading, confining SU {M{k + 1)) x SU{Mk) supersymmetric gauge 
theory. This theory approaches pure SU{M) SYM theory in the limit QgM <^ 1, while the dual 
supergravity description is reliable in the opposite limit, ggM ^ 1. 

In section 6 we connect the results of sections 2 - 5 to large Nc confining field theories such as 
pure YM. To leading order in l/A^, such theories are expected to reduce to free field theories of 
the gauge singlet bound states. The latter are expected to have a Hagedorn density of states at 
high mass, p{m) ~ exp(/3i^m). The entanglement entropy in such theories can be calculated 
by summing the contributions of the individual states. We show that this sum over states has 

^ This phenomenon has already been noted in one specific example - the static "AdS bubble," which is 
equivalent to the background of D3-brancs on a circle that we study in section 4 (we thank T. Takayanagi for 
pointing this out to us). 
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a very similar character to the finite temperature partition sum, with / playing the role of the 
inverse temperature f]. It converges for sufficiently large / and diverges below a critical value of 
/, since the large entropy overwhelms the exponential suppression of the contribution of a given 
state of large mass. In the thermodynamic case, this phenomenon is related to the appearance of 
a deconfinement transition. By analogy, it is natural to expect that here it signifies a transition 
between an entropy that goes like A'^^ at large I, and one that goes like A^^ below a critical value. 
Since the gravitational analysis reproduces this feature of the dynamics, we conclude that the 
system with QsN^ ^ 1 is in the same universality class as the one with QsNc ^ 1. 

In section 7 we comment on our results and discuss other systems which one can analyze 
using similar techniques. We also point out some general issues related to the proposal of pp. 



2 Holographic computation of entropy 

The gravitational backgrounds we will consider have the string frame metric 

ds^ = a{U) [/3{U)dU^ + dxf'dXf,] + gtjdy'dy^ (7) 

where (/i = 0,1,..., parametrize IR*^^^, U is the holographic radial coordinate, and 
(z = d + 2, ■ ■ ■ ,9) are the 8 — c? internal directions. The volume of the internal manifold, 

\{dy'^/d^g, (8) 

i=i 

and the dilaton, 0, are taken to depend only on U . 

The radial coordinate U ranges from a minimal value, t/o, to infinity. As f/ — f/o, a p-cycle 
in the internal (8 — (i)-dimensional space shrinks to zero size, so l^nt(f^) = 0. The vicinity of 
U = Uq looks locally like the origin of spherical coordinates in IR'''''^ (times a compact space), 
and we assume that all the supergravity fields are regular there. In particular, a{U) and ct){U) 
approach fixed finite values as f/ — > Uq. The fact that q;([/o) > implies that the string tension 
is non-vanishing. This is the gravitational manifestation of the fact that the dual gauge theory 
is confining. 

Examples of backgrounds in this class that will be discussed below are the geometries of 
coincident and Z)4-branes on a circle with twisted boundary conditions [10], in which the 
shrinking cycle is a circle {p = 1), and the KS geometry [H] in which it is a two-sphere {p = 2). 
In the D3-brane and KS cases, the dilaton is independent of U . 

We would like to use the proposal (E]) to calculate the entanglement entropy of A and B ([T]) 
in the geometry ([7]). Denoting the direction along which the line segment /; in ([1]) is oriented by 
X, the entropy per unit volume in the transverse IR*^"^ is given by 



dx^/H^^l + P{U){d.,UY (9) 



where we introduced the notation 

H{U) = e-'^V^d^ . (10) 



Due to the shrinking of the p-cycle, we have H{Uq) = 0. Thus, as U varies between Uo and 
oo, H{U) varies between and oo. It provides a natural parametrization of the radial direction 
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of the space ([7]). Near Uq, one has H ~ r^^, where r G [0, oo) is a natural radial coordinate, 
dr = ./p{U)dU. 

The quantity fllOp is simply related to the warp factor we get upon dimensionally reducing 
on the (8 — (i)-dimensional compact manifold. The resulting {d + 2)-dimensional Einstein frame 
metric may be written as 

dsl_^^ = k{U) [^{U)dU'^ + dxf'dXf,] . (11) 

A standard calculation shows that n{UY = H{U). It is a common assumption that the warp 
factor hi{U) is a monotonic function of the holographic radial coordinate. In particular, finiteness 
of the holographic central charge [12] , 

c~/?Ik^(k')"'' , (12) 

requires k to be monotonic. Since it goes to zero as f/ ^ f/o and to infinity as f/ ^ oo, it must 
be that > for all U . This implies H'{U) > 0, a fact that will be useful below. 

We need to find the shape U{x) that minimizes the action ([9]) subject to the constraint 
t/(x — >■ ±|) — > oo. Denoting by U* the minimal value of U along this curvejl and using the fact 
that the action does not depend directly on x, its equation of motion can be integrated once and 
written in the form 



a,f/ = ±^./j^M_i . (13) 



Integrating once more we find 



1{U*) = 2,/Hm r ^£^^= . (14) 
^ ^ ^ Ju* ^H{U) - H{U*) 



Plugging (fT3|) into ([9]) we find 

Sa _ 1 [^°° dU^/(3{U)H{U) 



Vd-i 2G!i°^yf/* ^HiU) - H{U*) 



(15) 



In the examples we study below, and probably much more generally, the integral in flMj) turns 
out to be convergent, while that in ( fTSi) is not. This is the reason for the appearance of the UV 
cutoff Uoo in the latter and its absence in the former. 

As mentioned earlier, the entropy Sa depends on the cutoff only via an / independent constant, 
which cancels in differences of entropies. This can be seen from (fTS!) as follows. Denoting by 
and U2 the solutions of (fT4l) for / = li and / = Z2, respectively, we have 



/oo 
dUsJ(3{U)H{U) 



H{U) J V H{U) 



(16) 



where we omitted an overall multiplicative constant and focused on the behavior of the integral 
in the UV region U 00. In that region H{U) 00, and we can approximate the integrand in 
m by 



SAih)-SAil2) - {H{Ul)-H{U;)) / dU^I§^^. (17) 



•^If the curve is smooth, this value is attained at x = 0, where d^U = 0. 
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The integrand in (fT7I) behaves as t/ ^ oo in the same way as that in (fT4|) . Thus, if the latter is 
finite and does not require introduction of a UV cutoff, the same is true for the former. 

To find the dependence of the entropy on / we need to determine U*{1) by solving (IHl), and 
then use it in ( ITSl) . In the next sections we will study specific examples of this procedure; here 
we would like to make some general comments on it. 

Consider first the limit U* — >■ oo. Physically, one would expect 1{U*) to go to zero in this 
limit since as / — > the minimal action surface should be located at larger and larger U. In 
terms of f|T^ this means that although the prefactor ^yH{U*) goes to infinity, the integral goes 
to zero faster, such that the product of the two goes to zero as well. We will see below that this 
is indeed what happens in all the examples we will consider. 

It turns out that / ( IT4|) also goes to zero in the opposite limit U* Uq. The prefactor 
-\/^ff(Z7*y goes to zero in this limit, and as long as the integral does not diverge rapidly enough 
to overwhelm it, / 0. Since any divergence of the integral as U* Uq must come from the 
region f/ ~ f/* ~ f/o, it is enough to estimate the contribution to it from this region. In terms 
of the coordinate r defined above, one has 



For p > 1 one finds that for small r^,, /(r*) ~ r*; for p = 1, /(r*) Inr*. In both cases, / 
in the limit — > (or, equivalently, U* Uq). 

We see that for small / the equation of motion 0131) has two independent solutions, one with 
large U* and the other with U* ^ Uq. The former is a local minimum of the action f|T5l) while 
the latter is a saddle point. We can interpolate between them with a sequence of curves which 
differ in the minimal value of f/, such that the solution with large U* is a local minimum along 
this sequence, while the one with U* ~ [/q is a local maximum. 

This implies that there must be another local minimum of the effective action, with U* smaller 
than that of the saddle point. This solution cannot correspond to a smooth f/(x), since then 
it would be captured by the above analysis. Therefore, it must correspond to a disconnected 
solution, which formally has U* = Uq, but is better described as two disconnected surfaces that 
are extended in all spatial directions except for x, and are located at x = ±|. 

The entropy corresponding to this solution is given by (see (IT^ ) 



By the above analysis it must be smaller than that of the connected solution with U* — Uq, but 
may be larger or smaller than that of the connected local minimum with large U*. 

We saw before that 1{U*), ( fT4l) . goes to zero both at large U* and as U* — ^ Uq. If the 
supergravity background is regular, one can show that between these two extremes / is a smooth 
function of U*, that remains finite everywhere. The simplest behavior it can have is to increase 
up to some point where 81/ dU* = 0, and then decline back to zero as U* oo. We will see that 
this is indeed what happens in all the examples we study below. 

Denoting the value of 1{U*) at the maximum by /max, this behavior implies that smooth 
solutions to the equation of motion (IT^ only exist for / < /max- As / — /max from below, the 
local minimum and saddle point discussed above approach each other, merge and annihilate for 




(18) 




(19) 



/ > /, 



max- 
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At first sight, the fact that there are no solutions to (fT3|) for / > /^ax n^ay seem puzzhng, but 
it is important to remember that this analysis only captures smooth connected solutions. As 
discussed above, for all I we have in addition a disconnected solution for which U'{x) is infinite. 
For / > /max the entanglement entropy Sa is governed by this solution and is given by ( |T9l) . For 
i < /max one needs to compare the entropies of the connected and disconnected solutions and 
find the smaller one. This difference can be written as 

2G|^ I^^Cconn) _ ^(disconn) j ^ J ^jj ^ I "^^^^^ ~^]~ j ^^V^ • (20) 



It is physically clear and easy to see from ( 1201) that for small / the connected solution with large 
U* has the lower entropy. As / increases, there are in general two possibilities. The connected 
solution can remain the lower action one until / = /max? or there could be a critical value /crit < ^max 
above which the right hand side of fl20|) is positive, so that the disconnected solution becomes 
the dominant one. In the first case there would be a phase transition at / = /max; in the second, 
the transition would occur at /crit, and in the range /crit < / < /max, the connected solution would 
be a metastable local minimum. In all the examples we study below it is the second possibility, 
^crit < ^max5 that is realized: as we increase /, the transition occurs before the connected solution 
ceases to exist. This is similar to the first-order finite temperature deconfinement transitions 
found in gravitational duals of confining gauge theories [10], [131 El [15] . 



3 D4-branes on a circle 



The low energy dynamics of D4-branes in type IIA string theory is governed by 4 + 1 di- 
mensional supersymmetric Yang-Mills theory with gauge group U{Nc) and 't Hooft coupling 
A = gsNJs- In order to reduce to 3 -|- 1 dimensions and break supersymmetry, we compactify 
one of the directions along the branes, x'^, on a circle of radius R4, x'^ ~ -|- 271R4, with twisted 
boundary conditions for the fermions. 

The low energy dynamics of this system, which was studied in [10] , depends on the dimen- 
sionless parameter A4 = X/Ra, and can be investigated using different tools in different regions 
of parameter space. For A4 <^ 1, it corresponds to pure Yang-Mills theory with gauge group 
U{Nc) and 't Hooft coupling A4 (at the scale -R4). In the opposite limit, A4 ^ 1, one can use a 
gravitational description in terms of the near-horizon geometry of the brane^ 



-24> 



-) 

R 



3/2 



R^/^u'/'dnl + (^^^ 



R^' 

u) 



3/2 



fiU)idx 



4\2 



(21) 



(22) 



where R is related to the five dimensional 't Hooft coupling via the relation R^ = ttX, and 



fiU) 



-'1 



4vr A 



(23) 



"'Here and below we set a' — 1. 
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As U ^ Uq, the radius of the x^-circle goes to zero; (?7, x^) form together a cigar geometry of 
the type described in the previous sections. 

Comparing (l2Til to (171) we identify a,P, Vint as, 

Vin. = ^(i?^f^)x2vri?4f^y^'v^=^^^^^V7M- (25) 
3 \RJ Qjji 



The combination (fTOl) is given in this case by 



(26) 

Note that H'{U) > for all U > Uq, as mentioned in the previous section. 

The explicit form of the background can be used to verify the assertions of section 2 about 
the behavior of 1{U*). In particular, it is easy to check that the integral f[T^ converges. For 
U* 3> Uo it is given by 

r ('^) 1 

l(U*) = 2i?3/2 X 2v^— >|4--= . (27) 

We see that / indeed goes to zero in the limit U* ^ oo, as expected. Similarly, one can check that 
it goes to zero in the opposite limit U* Uq. The full curve 1{U*) can be computed numerically 
and is plotted in figure 1. It has the qualitative structure anticipated in section 2. The maximum 
of the curve occurs at U* ^ 1.2f/o, with 

^ 1.418i?4 . (28) 

At larger values of /, there is no smooth solution to the equations of motion (fT3|) . 

Turning to the entanglement entropy Sa, following the discussion of section 2 we need to 
calculate the entropies of the connected solution fllSp and the disconnected one f ll9p . and compare 
them. The calculations of the individual entropies must be done with the UV cutoff Uod in place, 
but the difference of entropies is insensitive to it (see f|T6l) . fl20|) ). 

For the disconnected solution, the entropy can be calculated in closed form: 

- -^,-1/2^(10) l^oo - I^O) ■ (29) 

For the connected one it is given by f|T5l) . which in general has to be computed numerically. For 
small I one can again perform the integral using the fact that in this case U* ^ Uq. One finds 



(conn).. StT^ V,R^/^ . 

- q .-1/2^(10) I ^oo 

» Uq 



R^ 



(30) 



Comparing to (!29|l we see that for small / the connected solution has lower entropy, in agreement 
with the general discussion of section 2. The fact that the entropy (l30l) scales hke 1//^ at small 
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0^ 



Uol.2Uo 2 Uo 

Figure 1: 1{U*) for D4-branes on a circle. 



3U, 



I is indicative of 5 + 1 dimensional scale invariant dynamics. This is what one expects, since at 
short distances the dynamics on the wrapped D4-branes is described by the (2, 0) superconformal 
field theory in 5 + 1 dimensions. Indeed, we find 



S 



(conn) 



ri(disconn) 



-V2{2TlRi){2TlR 



32^/? 



10) 



— + 



(31) 



which is precisely the entanglement entropy of coincident M5-branes compactified on a circle 
of radius R^ and the M-theory circle of radius -Rio = Qs found in [1]. 

A naive use of (1291) . (!30l) suggests that the disconnected solution becomes the lower entropy 
one at / ~ R2/U0 ~ R4, not far from /max (1281) . Of course, the small I approximation leading to 
f l30|l is not valid there, and in order to determine the precise position of the transition we need 
to evaluate f|T5|) . The result of that evaluation is shown in figure 2, where we also exhibit the 
entropy of the disconnected solution and, for completeness, that of the saddle point discussed in 
section 2 as well. 

We see that, as expected, the saddle point entropy is larger than that of the connected and 
disconnected local minima for all /. It approaches that of the connected one as / — > /max, and the 
disconnected one as / — * 0. The entropies of the connected and disconnected solutions cross at 
I = Lit < ^max given by 



'■crit 



1.288i?4 



(32) 



As explained in section 2, the entropy is governed by the connected solution and exhibits non- 
trivial dependence on / for / < /crit, while for / > /crit it is governed by the disconnected one and 
is / independent (to leading order in 1/Nc). 
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4 D3-branes on a circle 

In this section we study the system of Nc D3-branes wrapped around a circle of radius -R3 with 
twisted boundary conditions for the fermions. The discussion is largely parallel to that of the 
previous section, and some of the results already appear in [7j, so we will be brief. 

Before the compactification, the low energy theory on the D3-branes is = 4 SYM with 't 
Hooft coupling A = QgNc. For finite R3 one finds at long distances a 2 + 1 dimensional confining 
theory. For A ^ 1 that theory is 2 + 1 dimensional YM with 't Hooft coupling A3 = A/-R3 [TU] . 
For A ^ 1 one can instead use a gravitational description in terms of the near-horizon geometry 
of the Nc D3-branes, 



h{U) 



u 



j] Ku){dx'f 



where 



= AttX 



u'- — 



(33) 
(34) 

(35) 



and the dilaton is constant, (piU) = 0. Comparing (133|1 to ([7]) we find 



U 



Uj h{U) 
The combination (fTOl) is given by 

H{U) = {2'K^R3fL^U%{U) . 



Mnt = 2Tl^R3L^U^/h{U) 



(36) 



(37) 
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It is again monotonically increasing with U, as expected. 

All the calculations of the previous section can be done in this case as well. The integral (HM 
is again convergent. For small / (and large U*) one finds 



(38) 



The extension to all U* is plotted in figure 3. The qualitative shape of 1{U*) is similar to the 
D4-brane case shown in figure 1. The maximum occurs at U* ~ 1.113f/o, and 



/max ^ 1.383i?3 . 

The entropy of the disconnected solution is given by 



I 



max 



2Uo 



Figure 3: /(t/*) for D3-branes on a circle. 



(39) 




^(disconn) ^ ^'RsLy ^Jj2^ _ f;2) 



2G 



(10) 
N 



(40) 



where Vi is the length of the strip. The entropy of the connected solution is exhibited in figure 
4. For small / one has 



^7°"^(/) 
Therefore, for small / we find 



(conn)^,^, _ TT^RsL'^Vi , 2 _a 

^00 



2G 



(10) 
N 



v^r(| 



2^^3 

3 



1/ 



(41) 



r((conn)/r\ Q(disconn) 



r(i, 



V,{27tR,)^ + ... 



(42) 
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^ (disconn) 




^cHt ^1 



max 



Figure 4: Entropies of the connected (blue and red) and disconnected (black) solutions for the 
wrapped D3-brane geometry. 



which is the entanglement entropy of the 3+1 dimensional A/" = 4 SYM theory compactified on 
a circle of radius [1] . 

As is clear from figure 4, the transition between the connected and disconnected solutions 
happens again at a value of / smaller than Zmax- The numerical evaluation gives 



/„it ~ 1.2376i?3 



(43) 



5 Cascading Confining Gauge Theory 



The background dual to the cascading SU{M{k + 1)) x SU{Mk) supersymmetric gauge theory 
is the deformed conifold X]^=i ~ warped by M units of RR 3-form flux. The relevant metric 



h~^'^{T)dx^'dx^ + h^'^{T)dsl , 



where ds^ is the metric of the deformed conifold 



ds. 



^^(rfr^ + (,T)+cosh^ (0 [igr + iar] 



Here 



and the warp factor is given by 



K(t) 



+ sinh 



(sinh(2r) - 2r) 



r 



1/3 



21/3 sinh 



r 



h{r) = {gMa'f2^'h-^'^ H dx ^'^'^'^^l \ sinh(2x) - 2x)^/^ 

sinh X 



(44) 



(45) 



(46) 



(47) 
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The dilaton is constant and we set it to zero. For the details of the angular forms Qi, see [TTl [T6] . 

The cascading gauge theory has a continuous parameter, gsM. The theory approaches the 
pure SU{M) SYM theory in the limit QsM — > 0, while the dual supergravity description is 
reliable in the opposite limit, QgM — > oo. In this limit the geometry describes a gauge theory 
with two widely separated scales: the scale of glueball masses, 

"^'^ = ' ^^^^ 

and the scale of the string tension at the IR end of space (the tip of the cigar), a/T7 ^ 

The metric dB]), is of the form (^^| with 
Using J gi A 92 A 93 ^ 94 A 95 = 647?^, we get 

Att^ 

V^, = ^ h'/^e''/'Ksmh\T) . (50) 
V6 

Thus, all the general formulae of section 2 apply, with U replaced by the standard deformed 
conifold radial variable r. 
We find 

H{t) = e-^'^Vl,a^ = —e^^/^h{T)K\T) sinh^(r) . (51) 

3 

H can be seen to be monotonically increasing with r as noted in section [2] from general consid- 
erations. The general equation f|T^ with these identifications gives /(r*) for the KS background. 
As in the previous sections, the integral is convergent. For large r*, we can approximate /(r) 
using the asymptotic forms valid at large r, 

h{r) ^ 21/33 (^?,Ma')'e"'/' [r - e'^^l^ , K ^ 1^l^e-^l^ , (52) 

H{r) n'e\9sMar - ^) - 2-2/V2/^(^.Ma') ^^e'^/^ (53) 

This leads to the simplified expression, 

/(r*) = 2i/3e"2/3^^Ma' H ^""'^'^^^ . (54) 



The main contribution is from the region r ~ r*; shifting r — > + ?/ and keeping the lowest 
order term in y we conclude that for large r*, 

= ^^^^^^^^e-/3,.Ma'v^e--/^ (55) 

r(i/6) 
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As earlier, / goes to zero as r* ^ oo. One can also verify, as outlined in section [2], that as r* ^ 0, 
/ goes to zero again. The full curve, computed numerically, is presented in figure [51 We see that 
it shows the same qualitative behavior as the other cases (figures 1,3). The maximum occurs at 
T* ^ 2.1 with 

/max ~ 1.00 mgii^ball • (56) 



1(T.) 



^max 





Figure 5: /(r*) for the KS geometry. 



We now turn to the entanglement entropy 5*^. As earlier, we have to calculate and compare 
the entropies of the connected ( ITSl) and disconnected ( fT9l) surfaces. As discussed in section 2, 
each of these entropies must be computed with a UV cut-off in place, but the difference of the 
entropies is UV finite. The entropy of the disconnected solution is found to be 



S 



(disconn) 



Vo 



22/3l67r3a'2 \ 2 



2roo/3 _ 



21 
— ( 
8 



,2roo/3 



2.194 



(57) 



where the finite additive constant was computed numerically. For the connected solution, we 
first consider an analytic approximation valid for small /: 



SaIt* 



(r- l/4)3/2eW3^, 



■22/3i67r3a'2 J^, ^^(r - l/4)e2- - (r* - l/4)e2-* 

Approximating this integral as we did for /(r*), we find 
M2e4/3 



(58) 



3v/^r(2/3) , 

2r(i/6) ^ ^ 



/3 



(59) 
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Thus, for / < 1/mgiuebaii, 



^(conn) _ ^(disconn) ^ -y. J^^MsJ^^ ^ log^ (m^.^ebau/) + ■ ■ ■ ■ (60) 



"327r3/2r (lY 



In the cascading theory the effective number of colors is a logarithmic function of the distance 
scale (mini [IE]: 

= ^QsM^ log(mgi,ebaiiO + • • • • (61) 
/vr 



We see that the finite piece of the entropy is 

27V?r^|fiVi(/) 

8r(i)^/ 



-V, 3 J +••• ■ (62) 



For a 3 + 1 dimensional conformal gauge theory, the finite piece of the entanglement entropy is 
indeed of the form N^{V2/P). Following [1], we may use a minimal surface in AdS^ x T^^ to find 
the entanglement entropy in the dual SU{N) x SU{N) SCFT [20] H 

Hence, the result (l62l) we find for the cascading gauge theory is a reasonably modified form of 
the conformal behavior. The same distance-dependent effective number of colors was found in 
evaluation of correlation functions in the cascading theory [T8l [19] . 

Going beyond the small / limit, we present the result of the numerical evaluation of Sa in 
figure [5l which compares the connected, disconnected and saddle point entropies. As expected, 
the saddle point entropy is always the largest and approaches the disconnected solution for small 
/ and the connected solution as / —>■ /max- The connected solution has the lowest entropy for 
small / and is the dominant contribution in this regime. The point at which the connected and 
disconnected solutions cross is l^it < /max, which is found to be 

/e.it ^ 0.95 m-ii,,,u . (64) 

For / > /crit, the 0{N^) entropy is /-independent as explained in section [2l 



6 Comparison to field theory 

It is natural to ask whether the transition at finite / that we found in confining gravitational 
backgrounds also occurs in large asymptotically free gauge theories, such as pure YM or 
J\f = 1 SYM with gauge group SU{Nc). The location of such a transition would have to be 
around the QCD scale, Icrit-^QCD ~ 1- At such scales the theory is strongly coupled and it is 
difficult to evaluate the entanglement entropy Sa directly. 

To proceed one can use the fact that, at large Nc, confining gauge theories are expected to 
reduce to free field theories of glueballs, whose density of states grows like 

p(m) ~ m'^e'^""' (65) 

^The extra factor of 27/16 compared to the result for AdS^ x comes from the fact that vol(r^^) = 
lfvol(55). 
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at large mass m. The inverse Hagedorn temperature Ph is of order 1/Aqcd, and a is a constant. 
Both are difficult to calculate from ffist principles. Most of the states that contribute to fl65l) 
are unstable resonances whose width goes to zero as Nc ^ oo. More generally, all interactions 
between the glueballs go to zero in this limit. At finite the spectrum ([25D is effectively cut 
off at some large mass scale. 

We can use the above picture to calculate the entanglement entropy at large Nc, by summing 
the contributions of the glueballs. To avoid UV divergences, we will consider the quantity 

dSAjl) 

C = l^j- (66) 

which, as mentioned in the previous sections, does not depend on the UV cutoff. Consider, for 
example, a free scalar field of mass m. It is clear that the non-trivial dependence of fl^Bl) on I 
is via the combination ml. We will be interested in the region ml ^ 1, where C{ml) can be 
calculated as follows. In 1 + 1 dimensions, the large / form of C{ml) has been obtained in [21] : 
it is given by 

C,(ml) = —IU{2ml) ^ (67) 
4 8 

A four-dimensional free scalar field can be thought of as an infinite collection of two-dimensional 
ones, labeled by the momentum in the transverse IR^, k, with mass m{k) = ^/mP~+^. Summing 
over these momentum modes we find the 3 + 1 dimensional version of (Ej 



CM) = —4- / d'kCi(m(k)l) ~ ' ^ ^ e-'""' . (6J 
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We see that the contribution of a single scalar field to the entanglement entropy is exponentially 
suppressed at large massif This is similar to the exponential suppression of its contribution to 
the canonical partition sum at finite temperature, with the role of the inverse temperature (3 
played here by 21. 

For a theory with a Hagedorn spectrum (l65l) of bound states, the total entropy is obtained 
by summing over all states. 



The integral converges for / > /9/f/2 and diverges otherwise. This is the analog of the usual 
Hagedorn divergence of the canonical partition sum at the Hagedorn temperature. There, the 
physical picture is that for temperatures below some critical temperature, that is believed to be 
somewhat below the Hagedorn one (221 EH], the system is in the confining phase and the thermal 
free energy scales like N^. Above that temperature, the system is in a deconfined phase and the 
free energy scales like N^. 

Similarly, for the entanglement entropy in gauge theory we expect that for / above some lent 
that is somewhat larger than /5h/2 the entanglement entropy is of order A*"^ and is given by the 
convergent integral fl69|) . while for / < l^nt the entropy is of order N^, in agreemeent with the 
divergence of (1^ . 

The resulting picture is qualitatively similar to what we got from the gravity analysis in 
sections 2-5. Of course, as usual, the details are expected to differ because in the gravity 
regime the theory contains two widely separated scales. One is the scale of the lightest glueball 
masses, which goes like I/-R4 in the D4-brane analysis of section 3, like I/-R3 in that of section 
4, and like mgiuebaii ( HHl) in the KS geometry. The other is the scale of massive string excitations 



living near the tip of the cigar, y/Tg, which is parametrically higher than the glueball scale. Since 

the exponential density of states comes from these string modes, we expect to be of order 



The transition point lent in the gravity regime is instead determined by the inverse of the 

— 1/2 

lightest glueball mass, and is parametrically larger than the Hagedorn scale Tg . Thus, as we 
decrease /, the transition at / = /crit to entangelement entropy of order happens long before 
Ph, ^crit ^ Ph- For example, in the cascading theory lent/ Ph ~ VglM. 

In the asymptotically free field theory regime, there is a single scale Aqcd and everything 
happens around it. One can interpolate between the two regimes by tuning the 't Hooft coupling 
(e.g. making g^M small in the KS example). Our results suggest that no phase transition is 
encountered along such an interpolation - the two regimes are in the same universality class. 

The arguments presented above apply directly to large theories. It would be interesting 
to investigate whether the phase transition we found continues to exist at finite Nc, and to 
characterize its order. Studying the entanglement entropy in pure glue SU(Nc) lattice gauge 
theory would therefore be very interesting. 

7 Discussion 

In this paper we applied the holographic method for calculating the entanglement entropy, intro- 
duced in [T] , to confining theories with gravity duals. In the simple case of entanglement between 

^The same is true for fermions and other higher spin fields. 




(69) 




1-1/2 



s 
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a strip of width / and its complement, we found an interesting phase transition as a function of 
/: for / < /crit the entropy is dominated by the action of a connected surface, while for / > lent by 
that of a disconnected one. After a subtraction of an /-independent UV divergent contribution, 
we conclude that the entropy is 0{N^) for / < /crit and 0(1) for / > lent- This transition is 
qualitatively similar to the confinement /deconfinement transition at finite temperature. 

Studying the thermal phase transition in confining gravitational backgrounds requires finding 
a SUGRA solution with an event horizon, and comparing its action with that of another solution 
which is horizon-free but has the Euclidean time periodically identified [10]. In general, these 
calculations are complicated and require a considerable amount of numerical work (see, for 
example, [131 [HI [U] ) • Studying the qualitatively similar transition for the entanglement entropy 
is much simpler; instead of finding new SUGRA solutions, one needs to find locally stable surfaces 
in previously known backgrounds. 

We also argued that a transition similar to the one we observed using the methods of [1] 
should occur in any confining large Nc gauge theory. This reasoning, and the several examples 
we have presented, make it plausible that any consistent gravity dual of a confining theory has 
to exhibit this phase transition. This is a useful prediction for any confining gauge/gravity dual 
pairs that remain to be discovered. 

The existence of the transition in the cases we have discussed is linked to p-cycles of the 
internal geometry that shrink in the IR. One could ask if this is the most general situation that 
results in the phase transition. As we showed, the monotonic function H{U) = e~^'^V^^a'^ has to 
vanish at the IR "end of space," U = Uq. On the other hand, a{Uo) should be non-vanishing for 
the string to retain its tension in the IR. This seems to restrict us to the vanishing of e~^'^V^nt- 
Thus, we should consider models where there are shrinking cycles and/or diverges in the IR. 

Curiously, one of the most widely used gravitational models of confinement [23], AdS^ with a 
hard IR wall at f/ = Uq, exhibits neither of these phenomena because both and \4nt are assumed 
to be constant. Therefore, for such a model the transition of the entanglement entropy does not 
seem to occur. This is not surprising, since the notion of the disconnected solution wrapping the 
entire geometry is not a priori well-defined in this case. A related problem is that the equations 
of motion are not satisfied at U = Uq, hence the boundary conditions are ambiguous there. 

There may exist a definition of the boundary conditions that allows the disconnected solution 
and produces a phase transition of the entanglement entropy (an encouraging sign is that the 
thermal deconfining transition does take place in the hard- wall model [25]). Indeed, when the 
hard wall model was considered in [1] the contribution from the part of the minimal surface 
lying along the hard IR wall was not included in the calculation; hence, it was treated as a 
disconnected surface. Justifying such a prescription may be a good problem for the future. 

Another popular phenomenological model is the "soft wall" model where space-time has the 
geometry of AdS^, while (piU) blows up in the IR [26] : 




Here, there is no shrinking internal cycle but the blow-up of the dilaton causes H{U) to rapidly 
approach zero at f/ = 0|!| In general, if H{U) ~ U^e-^/^^ as U ^ 0, one finds a finite /max 
(above which the connected solution does not exist) provided l3{U) has a pole of order 2g + 2 or 
less at U = 0. One can show this by similar means to those employed in section [2] where only 

^ For the soft-wall model a{U) = C/^, hence the string loses its tension at U = 0. However, the model is 
typically treated as a five dimensional field theory, so it is not clear if the string tension requirement needs to be 
imposed. 




dsl = U^ {U~^dU^ + dx^dx^) 



0(f/) = u- 



-2 



(70) 
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shrinking cycles were considered. In all the models considered in this paper so far, q = and (3 
had a pole of order less than 2. On the other hand, the soft wall model corresponds to g = 2 
while (3{U) = and hence still satisfies the criterion for the existence of a finite /max- For 

the soft-wall model one finds that there is indeed a transition between the disconnected solution 
stretching from U = to f/ = oo and the connected one that becomes unstable for / > /crit and 
stops existing at /max- 

We see that the entanglement entropy may be useful as a simple test of holographic models 
of confinement. More amibitiously, it would be nice to show that, if the confining background 
satisfies the supergravity equations of motion (neither the hard- wall nor the soft-wall do), then 
there is a phase transition of the entanglement entropy. 

Finally, it is important to understand the underlying reasons for the success of the geometric 
method of [1]. This prescription is designed to capture only the leading, 0{N^), term in the 
entanglement entropy. While it has a superficial similarity to probe brane calculations, it does 
not seem to be consistent to think of the bulk surface that appears in the construction as a 
brane. Indeed a brane with the worldvolume action iQ would have tension proportional to l/g^, 
and would back-react on the geometry at leading order in Qg. In any case, branes with the right 
properties do not seem to exist (see e.g. [27]). We need to formulate the problem in semiclassical 
gravity whose solution to leading order in G^"'* is the minimization problem proposed in [1]. 
Hopefully, this can pave the way to finding the 0{N^) corrections to the entanglement entropy 
and comparing them with field theory. 
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